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A scanning electron microscope (SEM) was used to measure the electron 
(minority carrier) diffusion length (LJ and the edge surface-recombi- 
nation velocity (V.) in zinc-doped Czochralski-grown InP wafers. Elec- 
tron-beam-induced current (EBIC) profiles were obtained in specimens 
containing a Schottky barrier perpendicular to the scanned (edge) 
surface. An independent technique was used to measure V,, and these 
values were used in a theoretical expression (Donolato, Ref. 12) for 
normalized EBIC. A fit of the experimental data with this expression 
enabled us to determine L a . 


I. INTRODUCTION 


The minority carrier diffusion length (L) is an important parame- 
ter in determining the performance of minority carrier devices, such as 
solar cells. In the past many different techniques have been used to 
determine L. Some rely on the measurement of the minority lifetime by 
means of photoluminescence [1,2]. The scanning electron microscope 
(SEM) is widely used for the measurement of L. With the SEM a high 
energy electron beam can be used in a line scan mode to generate a 
volume of charge carriers within the sample. The advantage of using an 
electron beam as opposed to other sources of excitation, such as 
optical [3] , is that the volume and depth of generation can be accu- 
rately controlled by varying the beam voltage. 

In the "normal collector geometry" the p-n junction or Schottky 
barrier is viewed edge-on. With the SEM in a line scan mode, the 
electron beam scans the semiconductor perpendicular to the potential 
barrier (Figure 1) . The generated charge carriers can then diffuse to 
the potential barrier where the electrons and holes are separated and a 
current, I(x,z), is generated in the external circuit. This current, 
referred to as electron-beam-induced current (EBIC) , reflects the 
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amount of excess carriers generated. The surface on which the beam 
impinges acts as a surface recombination path into which the generated 
minority carriers diffuse and are annihilated. We will refer to the 
recombination velocity of this surface as the edge-surface recombina- 
tion velocity, V, (cm/sec) . 

If we make the assumption that the electric field outside the 
junction space charge region is negligible, the transport of the 
generated minority carriers is purely diffusive. Under this condition, 
and if V,=0, the EBIC will decrease exponentially with increasing 
distance, x, from the junction as follows: 


I(x,z) - I(o)exp(-^) | 

L z-const. 


(1) 


(For a definition of the coordinate system refer to Figure 1.) A plot 
of log[I(x,z)] versus x, therefore, would result in a straight line 
from which L can be found. If the diffusion length is different on 
either side of the junction, the slopes will also be different. The 
plot will go through a peak which occurs at the metallurgical junction, 
as shown in Figure 2 [4] . In practice, however, V. cannot be neglected 
and the plot of log [I(x,z)] versus x is no longer linear but appears 
concave upward near the junction, becoming steeper with increasing V, . 

Increasing the electron beam accelerating voltage, V c , 
increases the depth at which the carriers are generated, thereby 
minimizing the effect of V,. However, for large beam voltages the 
electron range, R, may become comparable to the value of L, diminishing 
the resolution of the technique [5] . Figure 3 shows a plot of R as a 
function of the beam voltage, V 0 . Throughout this work V„ was limited 
to 15 KV. 

Several theoretical expressions have been derived for the induced 
current profile which incorporate the effects of V, . The first such 
expression was derived by Van Roosbroeck [6] from the solution of the 
diffusion problem for a point source of minority carriers at a depth, 
z, in a semi-infinite specimen. Other papers have been published which 
refer to Van Roosbroeck [7] . Expressions based on more realistic 
generation schemes, such as the uniform sphere or the spherically 
symmetric Gaussian, have also been derived [8-12]. These expressions 
give the induced current profile in terms of integrals of a modified 
Bessel function. In limiting cases (V.-.0, V.-.00) these lead to analyti- 
cal expressions of I(x,z) which have exact solutions. However, for 
arbitrary V., these have to be solved numerically . 


Donolato derived a simpler expression for I (x, z) through the use 
of the Fourier transform method [13] . The simplification results from 
the consideration of a two-dimensional study of the diffusion problem. 
Hakimzadeh et al . [14] have measured the minority carrier diffusion 

length in GaAs solar cells by fitting experimentally obtained EBIC 
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profiles to Donolato's equation. In their study, an independent tech- 
nique was used for the measurement of V, as a function of x [15] . 

In this work we have applied the experimental technique described 
in [14] to measure electron diffusion length, L nj in Zn-doped InP 
materials. A review of the theoretical approach is given in Section II. 


II. REVIEW OF THE THEORY 


The expression derived by Donolato is shown in Equation 2. 


I(x, z) 
1(0, z) 


- ^!—2 {exp{ ~^T~ ) ‘ 0.57xexp(-^i_nz c ) 

7T J 2 2 

x—^—erfc[— ([i-— ) ]}xsin(iDf) dK 

y/2 o 2 ( 2 ) 


where : 


k = 


L 


(3) 


_i 

= ( k 2 +X 2 ) 2 


(4) 


X = 


R 


(5) 


s = 


D 


( 6 ) 


where R is the range of the generation volume in (im, K is the wave- 
number, D is the minority carrier diffusion coefficient in cm 2 /sec, and 
I (0,z) is the maximum EBIC collected at the junction in amperes. If we 
assume the generation volume to be a three-dimensional Gaussian, the 
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range, R, will be given by [16] : 


R = 


0.027 6 AV a - €n 
(Z c ' 8S9 p) 




(7) 


where V 0 is in KV, A is the atomic weight in g/mole, Z is the atomic 
number of the target, and p is the density of the semiconductor in 
g/cm 3 . It can be seen that the two unknowns in Equation 2 are L and s 

The technique developed by Watanabe et 3.1. [15] was used to 

measure s as a function of x. Their Equation is repeated here in our 
notation : 


s = 


Ys 

D 


d 


In I{x, z) 




( 8 ) 


EBIC profiles were obtained along the same line scan for a number of 
accelerating voltages, from 4 KV to 15 KV. For each point, x, s was _ 
obtained from the slope of the ln[I(x,z)] versus Z 0 plot using Equation 
8 A typical plot is shown in Figure 4. The s values obtained in this 
manner were stoxed in an ASCII data file. The advantage of working with 
s rather than V, is that D need not be known. This eliminates errors 
resulting from a calculated value of D. 

To obtain L n , EBIC profiles for a 15 KV accelerating voltage were 
used. A program was written in FORTRAN to perform the integration in 
Equation 2 by approximating the generation volume by a Gaussian. The 
accuracy of this program has been checked previously [14] . When running 
the FORTRAN program the user is prompted for V of 1(0, z), the names of 
the ASCII data files containing experimentally obtained I (x, z) and the 
s values. I(0,z) was calculated by extrapolating the experimental plots 
of log [I (x, z) ] versus x back to the junction, as described in [14]. 

The program calculates values of L n for different points, x, and 
outputs these to an ASCII data file. 


III. EXPERIMENTAL DETAILS 


The aim of this work was to measure L n in InP materials. For this 
reason Schottky barriers were formed to minimize the effect of process- 
ing which may result from junction formation. We used Czochralski-grown 
zinc-doped InP wafers with reported carrier densities of 2 x 10 cm 
and 1 x 10 18 cm -3 , purchased from Crystacomm. All wafers were of (100) 
surface orientation with an uncertainty of 2° off axis towards (110) . 
Ohmic contacts were evaporated on the back (unpolished) surface. Gold 
Schottky contacts (2000 A thick) were evaporated on the front 
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(polished) surface in the form of 2.5 mm-diameter dots. The specimens 
were cleaved to expose the rectifying junction, as shown in Figure 5. 
SEM sample holders held these specimens in place and made electrical 
contacts to the front and back of the specimens. EBIC profiles were 
obtained by scanning the electron— beam along this cleaved edge, and the 
analysis was carried out as described in Section II. 

We observed an unexpected effect which prevented us from measuring 
L n in the lower doped materials. This effect and other results will be 
discussed in Section IV. 


IV. RESULTS AND DISCUSSION 


We found that in the lower doped specimens EBIC profiles did not 
decay very much with increasing x implying that charge collection 
occurs with approximate unity efficiency up to hundreds of micrometers 
from the Schottky barrier. An example is shown in Figure 6. This effect 
has been observed before in Si [17] and is thought to be due to an 
inversion layer which forms as a result of the interaction of the 
electron beam with the native oxide at the surface. The charge intro- 
duced by the electron beam on this surface is annihilated by the nearby 
holes and results in a negatively charged layer close to the surface. 
This inversion layer created near the surface results in the collection 
of the beam-generated minority carriers along the entire length of the 
scan. 


In Si [17] it has been shown that after about 20 successive scans 
the charge collection efficiency at large distances from the junction 
is progressively reduced [17], and finally normal behavior is re- 
established. This was not observed here for InP, even after many more 
scans . 

To alleviate this problem, we attempted to remove the native oxide 
prior to the SEM characterization by etching the edge surface in a 
solution of 10% HF [18] . However, even in the vacuum chamber of the SEM 
(which is - 10‘ ! or 10" 6 Torr) we were unable to keep the oxide layer 
from building up long enough to make the necessary measurements. 

This surface inversion effect was not observed in the higher-doped 
specimens since they require much more charge at the surface to create 
an inversion layer. EBIC profiles in these higher— doped specimens 
decayed with increasing x, as expected. The analysis described in 
Section II was applied to measure L n and s as a function of x in these 
specimens. Typical results are shown in Figure 7. 

Figure 7 illustrates that L n increases sharply as the junction is 
approached. This is due to the fact that the lateral extension of the 
generation volume is £ * 2a “ R/2 [13], therefore Donolato's equation 
is only applicable for x h R/2. For x •< R/2 the effect seen is an 
artifact since Donolato's equation no longer holds in this region. At 
large distances, x, L n is seen to increase again. In these regions the 
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EBIC measured is very small and comparable in magnitude to the beam 
current. This results in a large error in the measured values of L n . We 
therefore concentrate only in the middle regions where L n appears to be 
constant. In these regions the measured values of L n and s were as 
shown in Table I. These results are in the right ballpark as compared 
to previously published results which indicate that for p— type InP with 
a carrier concentration of 1 x 10 18 cm -3 L n is about 3 pm [19] . If we 
assume that D is a constant in these specimens and is approximately 
equal to 104 cm 2 /sec, the measured V, values ranged from 40.86 to 1.80 
x 10 3 cm/ sec . 

The aim of this work was to measure the L n values in the bulk. 
Therefore, it was necessary to ensure that the measured values were not 
in any way affected by surface defects. The edge surface of some 
typical samples were etched in concentrated HC1 acid for about 15 
seconds. This created a uniformly "pitted" surface as shown in Figure 
8. It has been shown that etching in concentrated HC1 for this length 
of time removes about 3 pm from this edge surface [20] . 

EBIC analysis was carried out on such treated samples, and L n and 
s values were determined as described in Sections II and III. Figure 9 
shows the results before and after treatment for a typical sample. It 
can be seen that although the s values have increased in the treated 
samples, L n appears to remain constant. 
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TABLE I . Measured values of L„ and s . 


Wafer Number 

Sample Number 

L n (^m) 

Range of s 
(xlO 3 cm' 1 ) 

3 

1 

0.64 

37.42 - 47.71 

3 

2 

1.36 

8.86 - 64.87 

3 

5 

0.78 

4.25 - 19.28 

3 

7 

0.39 

5.74 - 33.61 

3 

8 

0.45 

87.25 - 108.33 

3 

13 

0.44 

52.73 - 135.03 

3 

14 

0.77 

149.52 - 183.08 

2 

19 

0.26 

88.71 - 347.23 

2 

20 

0.38 

154.08 - 171.57 

2 

21 

0.36 

136.76 - 164.21 

2 

23 

0.26 

33.17 - 51.47 

2 

25 

0.39 

8.50 - 52.12 

2 

26 

0.73 

17.16 - 36.93 

2 

28 

0.34 

36.18 - 109 . 68 

2 

30 

0.50 

53.43 - 71.73 

2 

31 

0.39 

91.50 - 143.46 

2 

32 

0.50 

116.02 - 186.98 

2 

33 

0.78 

137.09 - 149.02 


All wafers had the following specifications: 

Crystacomm crystal #4420, zinc-doped with a carrier concentration of 
1 x 10 18 cm -3 . 
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Figure 7. Results obtained for four typical samples. 



Figure 8. Photograph of a pitted surface, created by etching the 
edge surface in HC1 . 
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DIFF. LENGTH. Ln (microns) £ OIFF. LENGTH. Ln (microns) 


(CRYSTACOMM 4420, SAMPLE #5) 



x (microns) 




